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Abstract. We prove optimal order error estimates for the Raviart-Thomas inter- 
polation of arbitrary order under the maximum angle condition for triangles and 
under two generalizations of this condition, namely, the so-called three dimensional 
maximum angle condition and the regular vertex property, for tetrahedra. 

Our techniques are different from those used in previous papers on the subject 
and the results obtained are more general in several aspects. First, intermediate 
regularity is allowed, that is, for the Raviart-Thomas interpolation of degree fc > 0, 
we prove error estimates of order j -I- 1 when the vector field being approximated has 
components in W^~^^'P, for triangles or tetrahedra, where < j < k and 1 < p < co. 
These results are new even in the two dimensional case. Indeed, the estimate was 
known only in the case j — k. On the other hand, in the three dimensional case, 
results under the maximum angle condition were known only for fc = 0. 
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1 Introduction 

The Raviart-Thomas finite element spaces were introduced in [20),l22j. and extended 
to the three-dimensional case by Nedelec [19] , to approximate second order elliptic 
problems in mixed form. After publication of that paper there has been an increas- 
ing interest in the analysis of these spaces and on the approximation properties 
of the associated Raviart-Thomas interpolation operator. This interest has been 
motivated by the fact that, apart from the original motivation, these spaces (or 
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rotated versions of them in the two dimensional case) arise in several interesting 
applications. For example in mixed methods for plates (see [3 [HI US]) and in the 
numerical approximation of fluid-structure interaction problems [S]. Also, it is well 
known that mixed methods are related to non-conforming methods [31 [TH] , there- 
fore, the Raviart-Thomas interpolation operator can be useful in some cases to 
analyze this kind of methods (see for example [T] where a non-conforming method 
for the Stokes problem is analyzed). 

The original error analysis developed in [^ \n\ is based on the so-called 
regularity assumption on the elements and therefore, the constants arising in the 
error estimates in those works depend on the ratio between outer and inner diam- 
eter of the elements. In this way narrow or anisotropic elements, which are very 
important in many applications, are excluded. 

For the standard Lagrange interpolation it is known, since the pioneering works 
[U [TBI EI] and many generalizations of them (see [2] and its references) , that the 
regularity assumption can be relaxed to a maximum angle condition in many cases. 

Error estimates for the Raviart-Thomas interpolation under conditions weaker 
than the regularity have been proved in several papers. In [1] the lowest order case 
k — was considered and optimal order error estimates were proved under the 
maximum angle condition for triangles and a suitable generalization of it for tetra- 
hedra, called regular vertex property. This result was extended in [15j to prismatic 
elements and functions from weighted Sobolev spaces. 

It is not straightforward to extend the arguments given in [Ij to higher order 
Raviart-Thomas approximations. In [11] it was proved that the maximum angle 
condition is also sufficient to obtain optimal error estimates for the case k = 1 
and n = 2 and in [M] that result was generalized to any fc > 0. Also in [14] . 
error estimates for any fc > and n — 3 were proved assuming the regular vertex 
property. 

The error estimates obtained in |14j require "maximum regularity" . To be 
precise let II^u be the Raviart-Thomas interpolation of degree fc of u on a triangle 
T then, it was proved in '14] that 

C 

||U - nfeU||L2(T) < /ly+Vli/'= + i(T) (1-1) 

sma ^ ' 

where we have used the standard notation for Sobolev seminorms, a and hx are 
the maximum angle and the diameter of T respectively, and the constant C is 
independent of T. However, an estimate like (jl.ip but with fc replaced by j < fc, 
only on the right hand side, cannot be proved by the arguments given in |14| and 
therefore a different approach is needed. Let us remark that this kind of estimates 
is important in many situations. In particular, the lowest order estimate 

C 

||u-nfcu||i2(T) < HtIuIhut) 

sma ^ ' 

is fundamental in the error analysis for the scalar variable in mixed approximations 
of second order elliptic problems. In particular the inf-sup condition can be obtained 
from this estimate (see for example [lO ] 112 ) ). 
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The maximum angle condition was originally introduced for triangles. For the 
three dimensional case two different generalizations have been given. One is the 
Kf I'zek maximum angle condition introduced in |17| : the angles between faces and 
the angles in the faces are bounded away from tt. Another possible extension is the 
regular vertex property introduced in [T] : a family of tetrahedral elements satisfies 
this condition if for each element there is at least one vertex such that the unit 
vectors in the direction of the edges sharing that vertex are "uniformly" linearly 
independent, in the sense that the volume determined by them is uniformly bounded 
away from zero. 

These two conditions are equivalent in two dimensions but not in three. Indeed, 
the Kfi'zek maximum angle condition allows for more general elements. This can 
be seen in the following way: consider the two families of elements given in Figure 
[U where hi, /12 and /13 are arbitrary positive numbers. Both families satisfy the 
Kf izek condition but the second family does not satisfy the regular vertex property. 




(a) (b) 
Figure 1 

Essentially these two families of elements give all possibilities. Indeed, the fam- 
ily of all elements satisfying the Kfizek condition with a constant -0 < tt (i.e., angles 
between faces and angles in the faces less than or equal to ^p) can be obtained trans- 
forming both families in the figure by "good" affine transformations (see Theorem 
12.21 for the precise meaning of this) . This result was obtained in [I] in the proof of 
Lemma 5.9. For the sake of clarity we will include this result as a theorem. On 
the other hand, the family of all elements satisfying the regular vertex property 
with a given constant (see Section 3 for the formal definition of this) is obtained 
by transforming in the same way only the first family in the figure. Therefore, to 
obtain general results under the Kfi'zek maximum angle condition (resp. regular 
vertex property) it is enough to prove error estimates for both families (resp. the 
Erst family) in Figure [T] with constants independent of the relations between hi, /12 
and h^. 

The error estimates in |14) for the general TZTk were obtained assuming the 
regular vertex property and the arguments given in that paper cannot be extended 
to treat the more general case of elements satisfying the Kfi'zek condition. On the 
other hand, as we have mentioned above, the arguments in I14j can not be applied 



3 



to obtain error estimates for functions in H^^^{T)'^ with j < k. For these reasons 
we need to introduce here a different approach. 

In this paper we complete the error analysis for the Raviart-Thomas interpo- 
lation of arbitrary order fc > 0. We develop the analysis in the general case of LP 
based norms, generalizing also in this aspect the results of previous papers. Our 
arguments are different to those used in previous papers. The main point is to 
prove sharp estimates in reference elements. 

Let us explain the idea in the two dimensional case. Consider the reference 
triangle T which has vertices at (0,0), (1,0) and (0,1). A stability estimate on 
T can be used to obtain the stability in a general triangle by using the Piola 
transform. Afterwards, error estimates can be proved combining stability with 
polynomial approximation results. 

The original proof given in [20] uses that 

||nfeu||^2(T) < c||u||^i(j^). 

In this way, the constant arising in the estimate for a general element depends 
on the minimum angle and so the regularity assumption is needed. The reason 
of that dependence is that the complete TJ^-norm appears on the right hand side. 
Therefore, to improve this result one may try to obtain sharper estimates on T for 
each component of n^u. Denote with Uj and Hkju, j — 1, 2, the components of u 
and its Raviart-Thomas interpolation respectively and consider for example j = 1. 
Ideally, we would like to have the estimate 

l|nfe,iu||^2(^) < c||wi||^i(^). 

However, an easy computation shows that if, for example, u = (0,a;|) then, IlfcU = 
^{xi,X2) and therefore the above estimate is not true. In other words, even for a 
rectangular triangle T, Ilfc^iu depends on both components of u. Now, the question 
is: which are the essential degrees of freedom defining Ilfc.iu? 

To answer this question one can try to "kill" degrees of freedom by modifying 
u without changing Ilfe.iu. A key observation is that if r = (0, g{xi)) then H^^ir = 
(we will give the proof of this result for appropriate reference elements in the 
three dimensional case). Therefore, if v = {ui{xi, X2), U2{xi, X2) — "2(2^1, 0)) then, 
IIfe,iV = Ilfc iu. But the normal component of v on the edge £2 contained in the 
line {x2 — 0}, i.e. V2, vanishes, and so do all the degrees of freedom defining 
Ilfc associated with that edge. Moreover, if we now modify the second component 
defining w ~ (ui(xi, X2), U2{xi,X2) — M2(a;i, 0) — X2a), for some a € Vk-i, we still 
have that W2 vanishes on £2 and that H^^iw = Ilfc^iu (because we are modifying 
V by adding a vector field belonging to the Raviart-Thomas space of order k). 
But, as we will see, it is possible to choose a in such a way that the degrees 
of freedom corresponding to integrals over T also vanish. Of course, it will be 
necessary to estimate some norm of a. We will give the details of the proofs in the 
three dimensional case. It is easy to see that the same arguments can be used to 
complete the arguments explained above for the two dimensional case. 

The new contributions of this paper can be summarized as follows: 
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• We prove error estimates under the maximum angle condition with order j + 1 
if the approximated function is in W^^^'P{T)"', n — 2,3, where < j < fc 
and 1 < p < oo. 

• Under the regular vertex property we obtain estimates of anisotropic type 
also for general fc > and 1 < p < oo. We also show that this kind of 
estimates is not valid under the maximum angle condition. 

Let us finally mention that the interpolation error estimates of anisotropic type 
are necessary when one wishes to exploit the independent element sizes hi, /i2 and 
hs to treat edge singularities in elliptic problems or layers in singularly perturbed 
problems. The dilemma is that such estimates hold, as we show, only for tetrahedra 
with the regular vertex property but it seems to be impossible to fill space by using 
this type of elements only. An anisotropic triangular prism (pentahedron) can, 
for example, be subdivided into three tetrahedra, from which only two satisfy the 
regular vertex property while the third is of the type of the element at right hand 
side of Figure 1 . The only known way out so far is discussed in [T5] . These authors 
use pentahedral meshes or tetrahedral meshes which are obtained by a suitable 
subdivision of a pentahedral mesh. Pentahedra based on a regular triangular face 
satisfy the regular vertex property by construction. For the approximation on 
tetrahedral elements they use a composition of two interpolation operators in order 
to avoid the above mentioned insufficiency with the tetrahedra which do not satisfy 
the regular vertex property. This approach is restricted to prismatic domains so 
far. 

The rest of the paper is as follows. In Section 2 we introduce notation and 
give some preliminary results on the conditions on tetrahedra that we will work 
with. Then, we prove stability in LP(T)^ for the Raviart-Thomas interpolation 
of arbitrary degree for functions in W^'P{T)^. These stability results are proved 
in Sections 3 for elements satisfying the regular vertex property, and in Section 
4 for elements satisfying the maximum angle condition. The estimates obtained 
under both hypotheses are essentially different but the results are sharp. Indeed, 
in Section 5 we show that anisotropic type stability estimates can not be obtained 
for the larger class of elements satisfying the maximum angle condition. Finally, 
in Section 6, we derive the error estimates from the stability results and standard 
approximation arguments. 

2 Notation and Preliminary Results 

In this section we recall some known results involving geometric properties of certain 
degenerate tetrahedra. Most of these results were proved in [I] and . 

Given a general tetrahedron T C IR"^, po will denote an arbitrary vertex and, 
for 1 < i < 3, £i, with \\£i\\ = 1, will be the directions of the edges sharing 
Po and hi the lengths of those edges. In other words, T is the convex hull of 
{po} U {po + hi£i}i<i<3. 



5 



We will use the standard notation for Sobolev spaces W'^'^{il) of functions with 
all their derivatives up to the order k belonging to LP{fl), denoting by |j • ||vv'=.j'(n) 
the associated norm. The same notation will be used for the norm of vector fields 
u G W'''P{il)^. As it is usual, we use boldface fonts for vector fields. 

With Vk{T) we denote the set of polynomials of degree less than or equal k 
defined over T C IR'^. The Raviart-Thomas space of order k is defined as 

TZTk = Vk{Tf + (xi,X2,X3)7'fe(r), 

and for u G Vl^^'^'(T)'^ the Raviart-Thomas interpolation of order k is defined as 
XlfcU e TZTk such that 

/ HfeU-npfe - I n npk ypkeVk{F), face of T, (2.1) 

JF J F 

[Uku pk-i = [ u pk-i ypk-ieVk^iiTf. (2.2) 

JT Jt 

In the rest of the paper the letter C will denote a generic constant that may 
change from line to line. 

Now we introduce the different conditions on the elements that we will use. 
The first one, called "regular vertex property" was introduced in [1]. 

Definition 2.1. A tetrahedron T satisfies the "regular vertex property" with a 
constant c > (or shortly, RVP{c)) if T has a vertex po, such that if M is the 
matrix made up with ^i, 1 < i < 3, as columns then \ det Af| > c. 

One can easily check that a regular family of tetrahedra (with the usual defini- 
tion of regularity given for example in f9j) verifies the regular vertex property. On 
the other hand, simple examples like that at the left hand side of Figure[T]show that 
arbitrarily narrow elements are allowed in the class given by RVP{c) for a fixed c. 

Despite the presence of anisotropic elements the regular vertex property arises 
as a natural geometric condition if one looks for Raviart-Thomas interpolation error 
bounds. Indeed, looking at the vertex placed at po, one can see that the family of 
elements satisfying RVP{c), have three normal vectors (those normals to the faces 
sharing po) uniformly linearly independent (see p] for more details). A reasonable 
condition, since the moments of the normal components of vectors fields are used as 
degrees of freedom in the Raviart-Thomas interpolation. Strikingly, as was shown 
in [1] , the uniform independence of the normal components can be somehow relaxed 
and error estimates valid uniformly for a wider class of elements can be still obtained 
for Ho (and for 11^ as we will show). More precisely, we will prove error estimates 
under the maximum angle condition defined below, which was introduced by Krizek 
in [17] and is weaker than the RVP. 

Definition 2.2. A tetrahedron T satisfies the "maximum angle condition" with a 
constant tp < ir (or shortly MAC{ip)) if the maximum angle between faces and the 
maximum angle inside the faces are less than or equal to ip . 
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Let us mention that the estimates obtained under RVP are stronger than those 
vahd under MAC. Indeed, in the first case the estimates are of anisotropic type 
(roughly speaking, this means that the estimates are given in terms of sizes in 
different directions and their corresponding derivatives). On the other hand, we 
will show that this kind of estimates are not valid for the more general class of 
elements verifying the MAC condition. 

The definition of the maximum angle condition, is strongly geometric. In order 
to find an equivalent condition, more appropriate for our further computations, we 
introduce the following definitions. In what follows, will denote the canonical 
vectors. 

Definition 2.3. A tetrahedron T belongs to the family T\ if its vertices are at 0, 
hiBi, ^1262 and ^.363, where hi > are arbitrary lengths (see Figure la). 

Definition 2.4. A tetrahedron T belongs to the family T2 if its vertices are at 0, 
hiei + h2e2, ^.262 and h^e^, where hi > are arbitrary lengths (see Figure lb). 

Note that elements in J^2 satisfy MAC(f ) but they do not fulfiU RVP{c) for 
any c. 

Lemma 2.1. Let T be a tetrahedron satisfying MAC^ip). Then we have 

1- If OL < (3 < ^ are the angles of an arbitrary face of T , then 7 > ^ and 
/3,7G [(^-V5)/2,V^]. 

2. If po is an arbitrary vertex of T and X ^ V' ^ 0,1'^ the angles between faces 
passing through po, then > ^ O'nd t/jjCf) £ [(tt — ^)/2,^]. 

Proof See [m. □ 

For a matrix M e H'^^^, ||M|| will denote its infinity norm. The arguments 
used in the following theorem are essentially contained in the proof of Theorem 7 
in [T71 page 516]. We include some details here for the sake of clarity. 

Theorem 2.2. Let T be a tetrahedron satisfying MAC{ip). Then there exists an 
element T £ !FiU T2 that can be mapped onto T through an affine transformation 
F{x) = MSc + c with \\M\\, ||M~-'^|| < C where the constant C depends only on ip. 

Proof. Given a tetrahedron T we denote with pi, i = 0, 1, 2, 3, its vertices and use 
obvious notations for its faces and edges. Let P0P1P2 be an arbitrary face of T 
and P3 its opposite vertex. We can assume that the maximum angle 7 of the face 
P0P1P2 is at the vertex po. Then from Lemma \TJ\ we have 

^ . / . TT- . -1 

sm 7 > TO := mm < sm — - — , sm ip > . 

Let ti and t2 be unit vectors along the edges PoPi and PoP2- We can also assume 
that the angle ui between the faces P0P1P2 and P0P1P3 is not less than the angle 
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between P0P1P2 and P0P2P3 (otherwise we interchange the notation between the 
vertices pi and P2). Then, again from Lemma |2. II we have 

sinLu > m. 

Now consider the triangle P0P1P3 and choose k £ {0, 1} so that the angle ^ at the 
vertex p^ is not less than that at the vertex pi-fc. Using again Lemma [2.11 we 
obtain 

sin^ > m. 

We take now ta as the unit vector along PfcPa and define A/q as the matrix made 
up with ti,t2 and as its columns. Since the columns of Mq are unitary vectors 
we have ||A/o|| < 3. Also, the adjugate matrix of Mq has coefficients with absolute 
value bounded by 2 and therefore, ||AfQ"^|| < 6/|detMo|. Then, to obtain the 
desired bound for ||M^^|| it is enough to show that | detMo| is bounded by below 
by a constant which depends only on ip. 

Consider the parallelepiped generated by the vectors ti,t2 and ts. Let z be 
its height in the direction perpendicular to ti and t2 and y the height of the face 
generated by ti and ts in the direction perpendicular to ti. 

Since ||ti|| = 1 we have 

I det Mq I = z sin 7 = ?/ sin ui sin 7 = sin ^ sin lu sin 7 > mP' . 

as we wanted to prove. 

Obviously, the same properties are satisfied by the matrix Mi made up with 
t2, — ti and ta, as its columns. 

Now, define hi = |poPi|, ^2 — IP0P2I and /13 — |pfcP3|. If fc = take T £ Ti 
with vertices at 0, ft-iei, /i2e2 and /laes and if fc = 1 take T £ J-2 with vertices at 
0, hiGi + /i2e2, ^.262 and ^1363. Then it is easy to check that x i— > Af^x + p^ maps 
T onto T. □ 

As mentioned above, the regular vertex property is stronger than the maximum 
angle condition. Indeed, the following theorem shows that, under _RVT'(c), the 
reference family in the previous theorem can be restricted to Ti. 

Theorem 2.3. Let T he a tetrahedron satisfying RVP{c). Then, there exists an 
element T £ Ti that can he mapped onto T through an ajfine transformation 
F(x) = Afx + po with ||Af II, ||A^~-'^|| < C where the constant C depends only on c. 
Furthermore, ifhi,i = 1,2,3 are the lengths of the edges of T sharing the vertex 
PO; we can take T £ Ti such that, for i = 1,2,3, hi is the length in the direction 

Proof. Let po and ii be as in the definition of _RVP(c) and hi be the length of the 
edge of T with direction ii. 

Take M as the matrix made up with ii as its columns. Since | det(Af)| > c and 
^i are unitary vectors then it is easy to check that ||Af|| < C and ||A^~^|| < C with 
a constant C depending only on a lower bound of | det(Af)| and therefore on c. 

Then, if T is the tetrahedron of Ti with with lengths hi in the directions e^, 
the affine transformation i^(x) = Afx + po maps T onto T. □ 
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Remark 2.1. It is not difficult to see that the converses of Theorems \2.2\ and \2.3\ 
hold true. Namely, the family of elements obtained by transforming Ti U T2 (resp. 
Ti)) by affine mapsSc 1-^ Mx+c, where \\M\\,\\M~^\\ < C, satisfies MAC{-4)) (resp. 
RVP{c)) for some ijj (resp. c) which depends only on C. 

3 Stability under the regular vertex property 

The goal of this section is to prove the stabihty in for the Raviart-Thomas 
interpolation of arbitrary order of functions in W^'P{T)^, for families of elements 
satisfying the regular vertex property. Precisely, the main result of this section is 
the following theorem. 

Theorem 3.1. Let fc > and T be a tetrahedron satisfying RVP{c). If Po is 
the regular vertex, £i,i = 1,2,3 are unitary vectors with the directions of the edges 
sharing po, hi,i = 1, 2, 3, the lengths of these edges, and Ht the diameter ofT then, 
there exists a constant C depending only on k and c such that, for all u G W^'P(T)^ , 
1 < p < 00, 



The theorem will follow by Theorem 12.31 once we have proved error estimates 
for elements in the family JFi. 

First we will prove appropriate estimates in the reference element T defined as 
the tetrahedron with vertices at (0, 0, 0), (1, 0, 0), (0, 1, 0) and (0, 0, 1). This is the 
object of the next two lemmas. Afterwards, estimates for elements in will be 
obtained by scaling arguments. 

We denote with Fi the face of T normal to n^, with ni = (— 1,0, 0),n2 = 
(0, —1, 0), na = (0, 0, —1) and ~ :^(1, 1, 1)- We will use the same notation for a 

function of two variables than for its extension to T as a function independent of 
the other variable, for example, f{x2,x^) will denote a function defined on Fi as 
well as one defined in T (anyway, the meaning in each case will be clear from the 
context). In the same way, the same notation will be used to denote a polynomial 
Pk on a face and a polynomial in three variables such that its restriction to that face 
agrees with pk. For example, for pk S VkiF^) we will write Pfe(l — X2 — x^,X2,xz). 
In what follows lik.i^ denotes the i-ih component of HfeU. 




Lemma 3.2. Let f e LP{Fi), g e LP{F2), and h e £^(^3). // 



u(a;i,a;2,a;3) = (/(a;2, X3), 0, 0), v{xi,X2,X3) = {0,g{xi,X3),0) 



and 



^{X1,X2,X3) = {0,0,h{xi,X2)) 
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then, their Raviart- Thomas interpolations are of the same form, namely, there exist 
Qi <= VkiFi), i ~ 1,2, 3, such that 

IlfcU = (gi(a;2,a;3),0,0), IlfcV = {0,q2{xi,X3),0), 

and 

IlfcW = (0, 0, q3{xi,X2))- 

Proof. Let us prove for example the first equality, the other two are obviously 
analogous. Since divu = 0, we have that divH/jU = and therefore, from a well 
known property of the Raviart-Thomas interpolation (see for example [71 [T^]), it 
follows that ftfcU e VkiTf. 

On the other hand, using now (j2.ip for i = 2, 3, and that U2 = uj, = 0, we have 

and then, taking — Ilfc.iU, we conclude that n^^iul^ = for z = 2, 3. Therefore 
Ilfe.iU = xiri for some n S Vk-i{T) and so, using now ()2.2p and again that U2 — 
U3 — 0, we obtain that, for i = 2, 3, Ilfc^jU = in T as we wanted to show. 

Finally, since divIIfcU = 0, it follows that ^ g^-^^" = and so, 11^ lU is indepen- 
dent of Xi . □ 

Lemma 3.3. There exists a constant C depending only on k such that, for all 

u= K, 7/2,^3) e w^'P{ff, 

\\fik,M\Lp(f) ^ C* (l|w»lli4'i.p(f) + l|divu||^p(j;)) , i==l,2,3. (3.1) 
Proof. From the previous lemma we know that, if 

V = (mi,U2 - U2ixi,0,X3),U3 - U3{xi,X2,0)) 

then, llfc lU — llfc iv. 

Let a,/3 e Vk-iiT) be such that 

/ {v2-X2a)pk-i = Q and / (ws - a;3/9) = ypk-i e Vk-iiT). (3.2) 

Observe that those a and /3 exist. Indeed, it is easy to prove uniqueness (and 
therefore existence) of solution of the square linear systems of equations defining 
them. 

Define now w = [vi, V2~X2a, ^3 — 0:3/3). Since (0, X2a, x^P) £ TZTk{T) it follows 
that llfc iv = life iv^r and therefore 11^ lU = Hj. iw. 
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Taking into account that ^2!^, =0 and ws\p^ — and the equations (|3.2p . it 
foUows that IIw is determined by the equations 

/ ftfc^iwpfc-i = / wipk-i Vpk-i eVk-i{f) 
Jf Jf 

j^Uk,2^Pk-i = \fpk-i eVk-i{f) 
^nfc,3wpfe_i = Vpfe_i e7'fc-i(f) (3.3) 
/ ft/c.iwpfc = wipk Vpfe e 7'fe(Fi) 

[ nfc.2wpfe = ypkeVkih) 
JF2 

f rifc.swpfc = ypkeVkih) 

/ (ftfe^iw + nfe,2w + nA;.3w)pA; = / {wi + + W3) pk ^pkeVkCFi). 

JFa JFa 



'Fi JFi 

Now, for pk £ Vk{T), we have 



divwpfe = - / w • Vpfc + / w-npfc 



/ w • Vpfe + ^ {wi+ W2 + W3)pk + w ■ npk 

Jf V 3 JFi Jdf\Fi 



but, from the definition of w, we have 



- /" w • Vpfc = - [ wi 1^ and [ w • npfe = - / wip^, 

Jf Jf Cl^l JdT\Fi JFi 

therefore, for all pk G Vk{T), 

—7= [ {wi + W2 + W3)pk = /"divwpfe+ [ wi^^+ [ wipk- (3.4) 
v3 J_F4 Jf Jf Ca;i Jp^ 

But, 

div w = divv — div (0, X2a, x^P) = div u — div (0, X2a, X3P). 
So, using (|3.4p . (|3.3[) . and standard arguments, we obtain 



nfc.iu||^p(j5j - \\iik^ivf\\^p(^fj 



< 



C(l|wi|lvi/i,p(f) + l|divu||^p(j.j + ||div (0,a:2a,X3/3)||^p(f) 
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Then, to obtain (|3.ip for i = 1, it is enough to show that 

|ldiv(0,a;2a,a::3/3)||^p(^) < C(||mi||^i,p(^) + ||div uH^p^^^). (3.5) 

For pk e Vk{T) we have 

= / (0, V2 - X2a, V3 - X3(3) ■ Vpk 

Jf 



div (0, V2 - X2a, V3 - 2:3/3) Pk + [(^2 - X2a)n2 + (ws - xsfi)?!^] pk- 

Jdf 

Now we take 2:2, X3) = (1 - xi - 2:2 - xz)pk-i with pk-i G 7'fc-i(r). ^hen, 

since pk = on F4, {v2 — X2a)n2 = on dT \ F4 and [v^ — x^P)!!^ = on dT \ F4, 
it foUows that, in the last equation, the boundary integral vanishes. Then, 



(1 -a;i -X2 -2:3) div {0,V2 -X2a,V3 -X3P)pk-i = 0. 
That is, for all pk-i G Vk-iif), 

{l-xi-X2-X3)div{0,X2a,X3j3)pk-i ^ / (1 - xi - 2:2 - X3) div (0, ■i;2, W3)pfc-i. 



Therefore, taking pk~i — div (0, 2:2a, 2:3/?) and applying the Holder inequality we 
obtain 



(l-2;i-2;2-X3) |div (0,2:20;, 2:3/3)1 < C||div (0, W2, W3)||^p^^-| || div (0, 2:2a, 2:3/3) 

But, since all the norms on Vk-i{T) are equivalent we conclude that, 

||div(0,X2a,2:3/3)||^p(j.) < C||div (0, U2, W3)|lip(f )■ (3-6) 
Now, observe that div (0, V2, W3) = div (0, U2, M3) and 

||div(0,M2,W3)||^p(y) < C (^\\divu\\^^^f^ + 

then, follows from ([glB]) . 

Clearly, the estimates for nfc^2U and nfc^3U can be proved analogously. □ 

From the previous lemma and a change of variables we obtain estimates for 
elements in J-'i. 

The Raviart-Thomas operators on the elements T and T will be denoted by 11^ 
and Ilk respectively. Analogous notations will be used for variables and derivatives 
or differential operators on T and T whenever needed for clarity. 
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Proposition 3.4. Let T ^ Ti he the element with vertices at 0, hiei, /i2e2 and 
^1363, where hi > 0. There exists a constant C depending only on k such that, for 
u (ui, -1*2, U3) G W^'P{Tf and i = 1, 2, 3, 



|nfc,iu| 



< c Win. 



«IIlp(t) 



/lilldivul 



LP{T) 



LP(T) 



Proof. Let u G W^'P{TY defined via the Piola transform by 
1 



G(x) 



detS 



i?u(x), X = i?x, with B 



hi Q 

/l2 
/l3 



It is known that (see for example [El ED]), 



and 



nfcu(x) = 



div u(x) 



1 



detS 
1 



detS 



Bnfcu(x), 
divu(x). 



(3.7) 



(3. 



Consider for example i = 1 (the other cases are of course analogous). Using (|3.ip 
we have 



p 

Lv(f) 



hih2h3 

Khi 



Ilfc.lU 



hih2h3 f „ 



< G|ll"illL(f) 



p 

LP(f ) 

P , lldi^lll 

3 



LP{T) 



dui 



p 



LP{T) 



as we wanted to show. 

We are finally ready to prove the main theorem of this section. 



□ 



Proof of Theorem \3 . 1\ To simplify notation we assume po = 0. From Theorem l2.3l 
we know that, if T S J^i is the element with vertices at 0, ft-iei, ^1262 and ^363, 
there exists a matrix M such that the associated linear transformation maps T onto 
T. Moreover, Me; = i = 1, 2, 3. 

Now, given u e M^^'^'(r)^ we define u e W^'P{TY via the Piola transform, 
namely, 

u(x) = — ^- — A'f u(x), x = Afx . 
^ ' detM ^ 

Using Proposition [33] after the change of variables x 1-^ x we have 



dvL 



dx-i 



LP{T) 



/i?.||divuF ~ 

t" "tP(T) 
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where hf is the diameter of T and denotes the vector (f|j, ffj)* • But, 

|i=detMM-|;i, divu(x).-^di;;S(X), (3.9) 
and hf < \\M^^\\hT. Therefore we arrive at 



|n,u||^,(^) < c\\M\n\M-Y Mi^^^^ + J^h^; 



3 a 

ou 



LP(T) 



+ /4||divu||^p(j,j 



and recaUing that ||M||, ||Af ^|| < C with C depending only on c, we conclude the 
proof. □ 

4 Stability under the maximum angle condition 

In this section we prove a stability result weaker than that obtained in the previous 
section but which is valid for families of elements satisfying the maximum angle 
condition. 

The estimate obtained here, although uniform in the class of elements satisfying 
MAC(V'), is weaker than the estimate obtained in Theorem 13.11 under the stronger 
ByP{c) hypothesis. Indeed, in front of each derivative, it appears the diameter Ht 
instead of the length of the edge in the direction of the derivative. However, our 
result is optimal. In fact, we will show in the next section that estimates like those 
in Theorem 13.11 are not valid in general under the maximum angle condition. 

The main result of this section is the following theorem. 

Theorem 4.1. Let k > and T be a tetrahedron with diameter hx satisfying 
MAC{ip). There exists a constant C depending only on k and "0 such that, for all 
ue W^^P{Tf, 1 <P< oo, 

||nfcU||ip(T) < C (||u||lp(t) + /lT||Vu||ip(T)) . (4.1) 



The steps to prove this theorem are similar to those followed in Section 3. Now 
our reference element T is the tetrahedron with vertices at 0, ei + 82, 62 and 83. 
For ni = (l,0,0),n2 = ^(1,-I,0),n3 = (0,0,1) and = ^(0,1,1) we denote 

with Fi the face of T normal to and with F2 the projection of F2 onto the plane 
given by a;2 = 0. 

Lemma 4.2. Let f e LP{Fi), g e ^^(^2), and h e ^^(^3). // 

u(a;i,a;2,a;3) = (/(a;2, X3), 0, 0), v(a;i, a;2, X3) = (0, ^(xi, X3), 0), 

and 

v/{xi,X2,x:i) = (0,0, ft.(a;i,a;2)) 
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then, their Raviart- Thomas interpolations are of the same form, namely, there exist 
qi <= VkiFi), i — 1,3, and q2 E Vk{F2) such that 



UkU = ((ji(x2,a;3),0,0), IlfcV = (0, 92(2;!, xa), 0), 
IlfeW = (0,0, 93 (xi, 0:2)). 



and 



Proof. The proof is similar to that of Lemma 13.21 We will prove the first equality, 
the other two follow in an analogous way. 

First, we have that divIIfcU = and therefore IlfeU e VkiT)^. Then, proceeding 
exactly as in the proof of Lemma [??^ we obtain that 11^ 3U = in T. Analogously, 
using now (|2.1|1 for i = 4 we have {Ilk.2'^ + Ilk,3u)\p^ — 0, and so 

nfe,2U + nfc,3U ^ {1-X2~ xs jr 

for some r £ Vk-i{T). Consequently, using now (|2.2p and that M2 = ii3 = 0, we 
obtain nfc^2U + nfc_3U = in T. Then, since we already know that nfe^3U = 0, we 
conclude that nfc_2U = in T. 

Therefore, IlfcU = {q, 0, 0) for some q e VkiT) but, since divIIfcU — 0, it follows 
that Ilfc.i is independent oixi. □ 

Lemma 4.3. There exists a constant Ci depending only on k such that, for all 



||nfe,iu||^2(f) < Ci + |ldivu|l^p(j;)) 



|nfc,2u||^2(f) < Ci{\\u2\\wi.P 



|nfe,3u||^p(J^j < Ci (^||w3||^i,p 



dui 








dxi 


+ 

LP(f) 


0X3 


LP(f )/ 



iivu|lLp(f)) 



(4.2) 
(4.3) 
(4.4) 



In particular, for i = 1,2,3, 



linfe,,u||^2(f) < C2 



( 



\ 



-'^llvi'i, p(T) 



E 









LP(f) 1 



(4.5) 



for another constant C2 which depends only on k. 

Proof. Let v = (ui, U2 — U2{xi,xi, X3), U3 — ^3(2:1, X2, 0)) and a, /3 G 7^fe-i(T) such 
that 



{v2 - {xi - X2)a)pk-i = and j ^{v^ - X3P)pk-i = ^Pk-i e 'Pfe-i(T), 

(4.6) 
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and define 

w = (-ui,M2 - U2(a;i,a;i,a:3) - (xi - 0:2)0, U3 - U3(a;i, X2, 0) - X3P). 
Then, since (0, (xi — 2:2)0!, x^P) G TZTk, it follows from Lemma that 

Ilfc^lU = Ilfc^iW. 

Now, taking into account the definition of w and (|4.6p . we have that II^w is defined 

by 



/ nfc_iwpfc_i 
Jf 


if 


.1 Vpfc_i e 7'fc_i(f) 


/ nfc,2wpfe_i 

Jt 


- ypk 


_i e T'fc-iCf ) 




= Vpfe 


-1 e ) 


JFi 


= WiPk 
JFi 


Vpfc e T'fclFi) 


/ (iifc^iw- nfc,2w)pfc 


= WiPk 

JF2 


Vpfe e Vkih) 


/ nfc,3wpfc 


= Vpfc 


e Vkih) 


/ (nfc,2w + nfc_3w) pk 


= {W2 + 
JFi 


wzjPk ypkGVk{Fi) 


But, using again (|4.6p we have, for e Vk{T), 




J^div {0,W2,W3)pk = - 


1 (0, W2, W3) ■ 


ypk 



(4.7) 



(W2n2 + W3n3)pk + -7= / (^2 + W3)pk 
dT\F4 V2 JF4 



(4.8) 



Then, it follows from and (giH]) that 

l|ni,fcw||^p(^^ < C(^||u;i||^yi,p(j,) + ||div(0,W2,W3)||^p(j,) 

< + ||divu||^p(j;j + ||div(0, (xi - X2)a,X3P)\\^p^f^'^ 

Therefore, to conclude the proof of (|4.2p it is enough to show that 

||div (0, {xi — 2:2)0, X3P)\\j^p^f-^ <C{\\ ). (4.9) 
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But, for all pk G Vk{T), we have 

(0,u;2, W3) • Vpfc - / div(0, W2,u;3) + / (^2/12 + it;3"-3)p/c- 
Jf Jdf 

Now, taking pk — (1 — a;2 — a;3)pfc_i with G Vk-iiT) the boundary integral in 
the last equation vanishes, and therefore we obtain 



{I-X2- a;3)div (0, {xi - X2)a, x-i(i)pk-i = / (1 - X2- a;3)div (0, M2, ^3)^, 



Then, (|4.9p can be obtained with an argument like that used for (|3.5|) . Clearly, the 
proof of inequality (|4.4p is analogous to that of (|4.2p . 

Now, to prove (|4.3p . take v = {ui — mi(0, a;2, a^s), M2, M3 — M3(xi, X2, 0)), a, /3 g 
■Pfc^iCf) such that 

/ (wi - a;ia)pfc_i = and / {v^ - X3P)pk-i = Vpfc-i e Vk-i{T), 
Jf Jf 

and define 

w = (wi - xia, t;2, 1^3 - 2^3/?)- 
Using again Lemma l4?2l and that (sia, 0, X3/?) G T^T^ we obtain 

nfe^2U = nA:,2W. 

In this case, it follows from the definition of w that IlfeW is defined by 

/ nfe,iwpfc_i = \fpk-i e Vk-i{f) 
Jf 

/ nfc,2wpfc_i = w2Pk-i Vpfe_i e 7'fe-i(f ) 

Jf Jt 

/ nfe,3wpfe_i = "ipk-i^Vk-iif) (4.10) 
Jf 

[ Hfc.iwpfc = ypk^Vk[Fi) 

JFi 

/ (life, iW - life, 2w)pfc = / {Wi-W2)pk Vpfe e 7^/0(^2) 

JF2 JF2 

I nfe.3wpfe = "^PkeVkiFs) 

JFs 

/ (nfe,2w + nfe,3w)pfe = / {w2 + w3)pk Vpfe e 7^/0(^4). 

But, it is easy to check by integration by parts that, for all pk G 'Pfe(T), 

div (0,U'2,U)3)pfe = - / W2-^ j= / 'W2Pk + —/= / {■W2+W3)pk (4.11) 

Jf 0x2 V2J% V2Ji?4 
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and 



/ div(wi,W2,0)pfe = - [ W2^^ + ^ [ W2Pk + ^ [ iwi~W2)pk (4.12) 

Jf Jf OX2 V2 V2 



Now, it follows from (gjUD, and ((TT^ that 

||n2,feW||^p(^) < C (|lw2|l^yi.p(y) + ||div(0,W2,W3)|lip(f) 

and therefore, using the definition of w, we obtain 

9m 1 



Idiv (wi,W2,0)|Lp 



LP(T) 



< C 



^2|| vi/i,p(T) 



dxi 



LP(T) 



LP(T) 



Then, to conclude the proof of (|4.3p we have to estimate the last two terms in the 
above inequality. From the definition of ^3 we have, for all pk € Vk (T) , 

dw3. 



du3 




d{xia) 




d{x3P) 




dx3 


+ 

LP(f) 


dxi 


+ 

LP{f) 


8X3 


LP(f )/ 



= 



dpk 
W3^ — 
T 0X3 



dx. 



-Pk 



W3n3Pk, 



dT 



but, if we take Pk — — X2 — X3)pk-i with pk-i G Vk-i{T) the boundary integral 
in the last equation vanishes, and therefore 



[l - X2 - X3)— pk^l 



It 9x3 
from which we obtain 



(I- X2- X3jT: — Pk-l 

0x3 



dix3f3) 



dx3 



In a similar way we can prove 



< C 



LP{T) 



du3 



dx3 



LP{T) 



d{xia) 


< c 


dui 




dxi 


LP(f) 


dxi 


LP{f) 



and so (|4.3p is proved. 



□ 



Proceeding as in the previous section we obtain now estimates for elements 
in J^2- 

Proposition 4.4. Let T G J^2 be the element with vertices at 0, hiSi + h2e2, ^1262 
and ^1363, where hi > 0. There exists a constant C depending only on k such that, 
for u = (ui, U2, U3) e W^^v{ff and i = 1, 2, 3, 



\lp(T) - ^ 



\Lp(T) 



3 



dui 



dxj 



LP{T) 







dxj 


LP{T) } 



(4.13) 
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Proof. We proceed as in the proof of Proposition l3.4l Recall that now our reference 
element T is the tetrahedron with vertices at 0, 81+62, 62 and 63. Therefore, the 
same linear map given by B in ProDOsition l3.4l maps T in T. Then, if u g W^'P{T)^ 
is defined via the Piola transform we have 



and 
Using 

fifeiU 



nfeu(x) = 

divu(x) 



1 



dots 
1 



detS 

and changing variables we have 



snfeu(x), 

divu(x). 



(4.14) 
(4.15) 



p 

LP(f) 



hih2h3 



< c 



hih2h3 



/ 

\\Ut 

V 



c 



V 



'LP{T) 



P 

LP{f) 

IP 

LP(f ) 



E 

dui 



dui 



LP{T) j=i 



LP{T) 



E 

j- 
3 







dxj 


LP(T) 1 



duj 




dxj 


LP(T) j 



and therefore (|4.13p is proved. 



□ 



Remark 4.1. For i = 1 and i = 3 a better result can be obtained. Indeed, by the 
same arguments used in the previous proposition, but using now {4--'^ '^i^d, li4-4^ , 
we can prove the following estimates, 



llfeiU 



LP(T) 



< c 



ILP(T) 



E^. 



duj 
dxj 



■ /i2l|divu| 



LP{T) 



LP(T) 



Anyway, this clearly depends on the particular orientation of the element and so, 
it does not seem to be useful for general tetrahedra. 

We can now prove the main result of this section. 



Proof of Theorem \4-l\ From Theorem 12.21 we know that there exists T e JFi U J-2 
that can be mapped onto T through an affine transformation x i— s- Afx + c, with 
||A/||, < C for a constant C depending only on -0. To simplify notation 

assume that c = 0. 

If T S then, T satisfies the regular vertex property with a constant which 
depends only on tp and so (|4.ip follows immediately from Theorem 13. II Therefore, 
we may assume that T G J-2 and has vertices at 0, hiBi + /i2e2, /i2e2 and h^e^, 
where hi > 0. 
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Given u e W^'P{T)^ we use again the Piola transform and define u € W^'P{T)^ 
given by 

u(x) = — - — Mu(x), x = Mx. 
^ ' detM ^ 

Tlien, using that 
changing variables and using ()4.f p in T we obtain 

iin,ur^,(^) < c\\M\nM-Y (iiur^,(^) +/z?,||Minii^ur^,(^)) 

concluding the proof. □ 

5 Sharpness of the results 

In view of the results of the previous sections, it is natural to ask whether the 
estimate obtained under the maximum angle condition could be improved. The 
goal of this section is to show that this is not possible. 

Consider the element T £ T2 with vertices at 0, hiGi + /i2e2, /i2e2 and ft.3e3 
and with diameter hx- We are going to show that the inequality 



LP{T) - ^ \ ll"llLP(T) ~^ ^ "-3 



infc.25IU„,^. <c\ l|SL„,^. + > h 



^ diii 



hT\\diyu\\^p,f. , (5.1) 



LP(T) 

LP(T) 



with a constant C independent of hi,h2 and /13, does not hold for some u = 

(wi,M2,%) e w^'P{Tf. 

Suppose that (|5.ip . with C independent of /ii,/i2 and /13, holds true for all 
u G W^'P{T)^ . Let T be the reference element used in SectionUl i.e., T has vertices 
at 0, ei + 82, 62 and 63. Then, with u G W^'P{ff we associate u G W^-P{f)^ 
defined via the Piola transform with the linear transformation used in the proof of 
Theorem 14.11 

To simplify notation we drop the hat from now on and write u instead of u and 
Xi the variables in T. 

A simple computation shows that from inequality (|5.ip we obtain 



3 



\\Uk,2U\\^^f^f^ < |^_^/l,;||u,||,^i,^(j?) +/lT||divu||^p(f) 

Then, taking hi = = (with h2 < 1), we would have 

\\Iik,2u\\Lp^f) ^ ^ (^2||ui|lvi/l,p(f) + I|w2|lvyl,p(f) +^2||M3|lvyl,p(f) + II div u|| ^^^^^ 

and letting /i2 — > we would arrive at 

l|nfe,2u|l^p(f) < c(||w2|lvi/i,P(f) + l|divu||^p(j,)) . (5.2) 
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However, we are going to show that there exists u e W'^'P{T)^ for which inequahty 
(|5.2p is not vahd. In fact, in the next proposition we will give, for each fc > 0, a 
function u G W^'P{T)^ such that the right hand side of (|5.2p vanishes while the left 
hand side does not. We will use the notation of Section [3] for the faces of T. 

Proposition 5.1. For k > 0, the function u(a;i, 0:2, a;3 
verifies div u = 0, U2 = and llk.2U 7^ 0. 

Proof. We consider the case fc > 1 (the case fc = follows analogously). Since 
divu = we have Hfe^iu, Ilfc^au e Vk{T). 

Now, using that mi = on Fi and U3 = on Fg, it follows from the definition 
of IlfeU that 



and 



/ Uk,3npk = ypkeVkih). 



Then Ilfc^iu = xia and Ilfc.aU = x^P with a, (3 ^ Vk-iiT). Also from the definition 
of IlfcU we have 

/ (Ilfc^iU - nfc,2u)pfc = / (mi - U2)pk ypk G Vk{F2), 

and then, if nfe,2U — 0, we would obtain 

x^{(x-x\)pk = ^ VpkeVk{F2). 

F2 

But taking pk — a{xi, xi, x^) — x\ in the last equation, we arrive at a(a;i, a;i, 2:3) = 
x\, but this contradicts the fact that a e Vk-i{T). Then we have H;;. 2U 7^ 0. □ 



6 Error estimates for RT interpolation 

We end the paper giving optimal error estimates for Raviart-Thomas interpolation 
of any order. These estimates are derived from the stability results obtained in the 
previous sections combined with polynomial approximation results. 

Let us recall some well known properties of the averaged Taylor polynomial. 
For a convex domain D and any non-negative integer m, given / g M/f'™+^(D) the 
averaged Taylor polynomial is given by 

Q™/(x) = ^y^T„J(y,x)dy, 

where 

r„J(y,x)= ^ i?"/(y)iEi^Z)^. 

|a|<m 
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Then, there exists a constant C, depending only on m and D (see for example 
[51[T1]), such that 



WD'Hf- 



2/)I|lp(D) < C 



il+i2+i3=in+l 



Qm+1 j 



(6.1) 



LP{D) 



whenever < \[3\ < m + 1. 

As a consequence of these results we have the following approximation result for 
elements satisfying the regular vertex property. Given a function /, _D'"/ denotes 
the sum of the absolute values of all the derivatives of order m of /. 

Lemma 6.1. Let T be a tetrahedron satisfying RVP{c) such that po is the regular 
vertex, £i,i — 1,2,3 are unitary vectors with the directions of the edges sharing po, 
hi,i = 1,2,3, the lengths of these edges, and Ht the diameter of T . Then, given 
u G Vl/"+^'P(r)3, m > 0, there exists q e V,niTf such that. 



5(u - q) 



LP(T) 



gm+i 



u 



(6.2) 



LP(T) 



(6.3) 



and analogously for '^^"^ '^'^ with j ~ 2, 3. Also, 

lldiv (u - q)IUp(T) < C/j??|p™divu|Up(T) 
where the constant C depends only on m and c. 

Proof. To simplify notation we assume again that po = 0. From Theorem 12.31 we 
know that there exists a matrix M such that its associated linear transformation 
maps T onto T, where T is the element with vertices at 0, hiei, h2G2 and h^e^,. 
Moreover, the norms of M and of its inverse matrix are bounded by a constant 
which depends only on c. 

Now, as in the proof of Theorem [3711 we define u e M^™+1'P(T)^ via the Piola 
transform and 

QmU = {QmUl,QmU2,QmU3) G Vm{T)^ , 

where QmUj denotes the averaged Taylor polynomial of uj. 

Using the estimate (|6.1|) on the reference element T which has vertices at 0, ei, 
62 and 63, and a standard scaling argument we obtain 



d{u~ Q™u) 



dxi 



Lp(T) 11+22+13=™ 



Qm+1 



U 



dxl'+^dxl^dx'i 



LP{T) 



Then, defining q €E VmiT)^ via the Piola transform, that is. 



1 



-MQ™u(x), x = Afx, 



det M 

follows by changing variables as in the proof of Theorem 13.11 
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On the other hand, since 

divQmU Qm-idivu, 
using again (j6.ip in T and a scahng argument we obtain, 

|ldi^(S-Q,„5)||^,(f) < C/i^||5"di^u||^,(^) 



and therefore, (j6.3p foUows by using the properties of the Piola transform stated in 

(13:91). □ 



We can now state and prove optimal error estimates for elements satisfying the 
regular vertex property. Our theorem generalizes the results proved in [T], where 
the same error estimate was proved in the case fc = 0, as well as those proved in 
|14j . where the estimate was proved for any fc > but only in the case m = k. 

Theorem 6.2. Let k > and T be a tetrahedron satisfying RVP{c). If Po is the 

regular vertex, £i,i = 1,2,3 are unitary vectors with the directions of the edges 
sharing po, hi,i = 1,2,3, the lengths of these edges, and Ht the diameter of T 
then, there exists a constant C depending only on k and c such that, for < m < k, 
1 < p < oo, andue W^™+i'P(r)3, 



u 



< C 




d 



im+l. 



LP{T) 



de-^dt^dPi 

Proof. Since m < k, for any q € Vm{T)'^ we have 

u - rifcU = u - q - Ilfc (u - q) 
and therefore, applying Theorem 13. 11 we obtain 
llu-HfcU _ 

d{u. 



D"^divu\\Lp^T) 



< C 




It 



/iT II div (u — q) 



LP{T) 



\LP(T) 



Then, taking q e Vm{TY satisfying (|6.2p and (|6.3[) we conclude the proof. 



□ 



Also optimal error estimates under the maximum angle condition can be proved. 
We state the results in the following theorem. 

Theorem 6.3. Let k > and T be a tetrahedron with diameter hx satisfying 
MAC{ip). There exists a constant C depending only on k and "0 such that, for 
0<m<k,l<p<oo, and u e VF^+i'P^T)^, 



LP(T) 



<C/j™+^||n-+^u|Up(T) 



The proof is analogous to that of the previous theorem, using now the stability 
estimates obtained in Theorem 14. 11 and so we omit the details. 
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